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ABSTRACT 

We describe a complete volume limited sample of nearby active galaxies selected by their 14- 195 keV 
luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We 
describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy 
properties. The active sample appears to have no bias in terms of AGN type, the only difference 
being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. 

In the luminosity range spanned by the sample, logLi 4 _i 95 fc e y[e? , gs^ 1 ] = 42.4-43.7, the optically 
obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant 
spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does 
not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, 
comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the 
optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the 
X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert Is 
can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) 
emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges 
the gap to the dusty torus. 

Subject headings: galaxies: active - galaxies: Seyfert - galaxies: nuclei - X-rays: galaxies 


1. INTRODUCTION 

Feeding and feedback have become a paradigm of 
galaxy evolution models: by quenching star forma¬ 
tion, feedback from active galactic nuclei (AGN) is 
thought to shape the galaxy luminosity function and 
create the bi-modal colour sequence in galaxy p opula- 
tions (jBenson et, al.l l200Bt iKauffmann et, all 1200311 But 
the prescriptions used in models (iSpringel et al.l 120051 : 
iCroton et al.l 1200(1 iSomerville et all 1200811 are relatively 
simple because observations have focussed on the ques¬ 
tion of where the inflowing material originates (mergers 
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versus secular evolution) and on integrated galaxy prop¬ 
erties. To redress this, we need, in addition to large 
scale data, also to understand how the material flows to¬ 
wards the black hole (BH) on smaller scales. However 
these scales cannot be spatially reso l ved a t z > 1 where 
co-evolution largely occurs dFabianl [20121 1. Instead, it 
is local galaxies that currently offer the only opportu¬ 
nity to guide the small scale aspects of galaxy and BH 
co-evolution models. That many local AGN - in par¬ 
ticular Seyferts - are disky systems does not necessar¬ 
ily reduce their validity as templ ates for m ore distant 
galaxies. The nested simulations of lHopkins fe Quataertl 
(120101! provide an important insight in this respect: they 
suggest that as one looks further inside the central kilo- 
parsec, disk processes become increasingly important in 
driving gas inwards, independent of what has occurred 
on large scales. Observations also appear to confirm that 
even at z = 1-2 disk processes are a major contributor 
to AGN fueling feocevski et al.l 120121: [S chaw inski e t al 


2014). 


2012; IKarouzos. Jarvis, fe B onficld 2 014 IVillforth et al 


In nearby archetypal objects, integral field spec¬ 
troscopy (at optical and infrared wavelengths, and 
sometimes with adaptive optics) has probed kinemat¬ 
ics on 10 -100 pc scales, lead i ng to insights in both in¬ 
flow (see IStorchi -Bergmarml l2014aT) and outflow (see 
IStorchi-Bergmannl I2014bfh Despite the inherent com¬ 
plexities of performing detailed systematic studies with 
such techniques, and because of the large amount of 
observing time needed, there has been some progress 
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Miille r-Sanchez et al. 1201 It iRupke fc Veilleuxl I2011L 
20131 iMiiller- Sanchez et all I2015H . and in a few cases 


with a matched sample of inactive galaxies for compari- 


son 

Dumas et al. 2007!: IWestobv et al.i 2012: Hicks et al.l 

2013 

: Davies et al. 20141). There are also two larger stud- 


ies of local AGN based on integral field spectroscopy, 
but which do not address spa t ially r esolved structures 
or kinematics: iBurtscher et al.l (l2015l l focus on the spa¬ 
tially unresolved near-infrared non-stellar continuum in 
Seyferts and other low luminosity AGN (using the spatial 
information o nly t o hel p characterise this co mponent); 
iSturm et al.l (|2011h and iVeilleux et al.1 (120 1 .'11 1 focus on 
ULIRGs and QSOs, for which the data are spatially un¬ 
resolved due to the long wavelengths observed. 

In this paper, we present a complete, volume limited 
sample of nearby bright hard X-ray selected AGN, in a 
luminosity range that overlaps with AGN at higher red- 
shift. This statistically meaningful sample of 20 AGN is 
complemented by a sample of inactive galaxies that are 
matched in mass, morphology and inclination. In Sec¬ 
tion [2] we describe the scientific rationale for, and the 
selection of, both the active and inactive samples; and 
outline the observations that have been started in or¬ 
der to spatially and spectrally resolve the processes that 
drive and regulate gas inflow and outflow on small scales. 
We then look at the AGN properties in Section [3l in par¬ 
ticular the fraction of Sy Is and Sy 2s, and the neutral gas 
columns. Reconciling this with the numbers reported in 
the literature leads to a discussion of optical obscuration 
and X-ray absorption. Section [4] focusses on the issue of 
X-ray unabsorbed Sy2s. Then in Section [5] we look at 
what the differing luminosity dependencies of the opti¬ 
cally obscured and X-ray absorbed fractions can tell us 
about the torus and broad line region (BLR). We have 
chosen this approach - avoiding the need to invoke torus 
models - because current torus models are in n eed of 
additi onal basic constraints. As pointed out by Honig 
(l2013f) and lFeltre et al.1 (12012H . different models can come 
to contradictory conclusions and model parameters are 
often degenerate. 

2. SAMPLE 

Most of the studies in the literature that focus on the 
structure and processes occurring in the circumnuclear 
regions of samples of local AGN, have selected objects 
using the optical classification (e.g. as Seyfert or LINER) 
as the only defining characteristic of the AGN. An ob¬ 
vious reason for doing so is that there is a wide vari¬ 
ety of local AGN and such studies often try to span a 
range in AGN types. However, this may lead to confus¬ 
ing results because of the impact of accretion ra te, or 
equivalently and as a proxy for it, lumino sity (|Hol 120081 : 
IDavies et al .11 2 0121 : iKormendv fc Holl20 1 3f) . For example, 
a bright Seyfert such as NGC 1068 with log Lagn ~ 410 

a Throughout this paper we use units of erg s 1 for luminosities. 
Lagn refers to the bolometric luminosity of an AGN. L 14—195 
is the observed luminosity in the 14-195 keV band. 1 / 2—10 and 
^ 2—*10 are respectively the observed and absorption corrected (i.e. 
intrinsic) luminosities in the 2-10 keV band. Z/i 2 /xm refers to 
the monochromatic nuclear 12 /im luminosity, measured at sub - 
arcsecond spatial resolution (for details see lAsmus et al.l 120141) . 
We also adopt the relations log Lagn = 1-12 log Z/i 4—195 — 4.23 
(| Winter et al.l 120121 ) and log— 1-06 log L 14-195 — 3.08 
(|Winter et al.ir2009l) . 


requires an accretion rate of order O.lMQyr -1 , while a 
low luminosity AGN such as M81 with logL^Giv ~ 41 
will have an accretion rate ~ 10 -5 MQyr -1 that is 
4 orders of magnitude lower. This differ ence has a 
drama tic i mpact on fuelling requ irements. iHicks et al.1 
(|2009f ) and IMazzalav et al.1 (120131 ) show that one can ex- 
pect to find a gas mass of order 10 7 M 0 within a ra¬ 
dial scale of ~ 10 pc. Assum ing a 1% efficiency to bring 
this to accretion disk sc ales (|Miiller Sanchez et al.l[2009c 
iSchartmann et aEl I2010T) . this mass alone is sufficient to 
fuel M81 for lOGyr, but it would supply NGC 1068 
for only 1 Myr. In this context, an accretion rate sim¬ 
ilar to that for NGC 1068 is clearly very significant and 
may only be sustained via a relatively efficient inflow, 
perhaps driven by a coherent dynamical mechanism on 
scales of 0.1-1 kpc whe r e the gas reservoir is 1 0 8 -10 9 M Q 
(| Sakamoto et all 119991 : IMazzalav et al.112013T) . In con¬ 
trast, low accretion rates comparable to that for M81 
could be supplied either by a gas reservoir on small scales, 
or by rather inefficient inflow from circumnuclear scales. 
Observationally, such an effect is implicit in the relation 
bet ween sta r formation age and accretion rate presented 
by IDavies et al.1 (120071! : while AGN with higher lumi¬ 
nosities (~ 10 45 ergs -1 ) or Eddington ratios (0.1 1) may 
be associated with young post-starbursts, this is not the 
case for AGN with luminosities two orders of magnitude 
lower. It could al so help to explain the results of a survey 
of molecular gas ()Haan et al.ll2009t lGarcia-Burilldl201l( l 
which found that gravitational torques could drive gas 
inwards for not more than half of their sample. And 
the role of AGN luminosity has been highlighted in the 
context of molecular outflows, in the s ense th at onl y very 
luminous AGn drive massive outflows (|Sturm et al.ll201lb 
IVeilleux et, alJl2013h . 

Our premise for defining a new sample is that the black 
hole accretion rate, traced by the AGN luminosity, is a 
key issue that needs to be considered in order to provide 
a context both for the mechanisms driving and regulating 
inflow and outflow as well as the derived flow rates. To 
address this issue, our sample is based on bright local 
AGN in which the accretion rate is relatively high. We 
plan to analyse them in the context of matching inactive 
galaxies, for which the accretion rate is, by definition, 
orders of magnitude lower. 

An obvio us co n cern here is t he timescale _of AGN 
varia bility (jNovak. Ostriker. fo Giottll2011llHickox et al.1 
120141 1 with respect to the timescale of the phenomenon 
(e.g. starburst, dynamical process, etc) that is being as¬ 
sessed in the context of accretion onto the black hole. 
If an inactive sample is used naively as a control in a 
direct comparison to the active sample, the implicit as¬ 
sumption must be that any link to accretion is strong 
enough to persist even when averaged ove r the lifetime 
of the phenomenon being studied. However. IDavies et all 
(I2014H argue and show that with a more judicious use of 
an inactive sample, one can make progress even if this 
assumption is not strictly met. 

Buil ding on the smalle r sample of IHicks et al.1 (|2013f l 
and IDavies e t, al.l (I2 014H. and also o n the work sum¬ 
marised bv iStorchi-Bergmannl d2014al ibil. we plan to ad¬ 
dress a few key questions in a statistically robust way: 

1) Does star formation play a decisive role in either 
driving or hindering gas inflow to AGN? This question 
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Fig. 1.— Plot of the redshift versus luminosity for AGN in the 58-month Swift-BAT survey that are readily observable from the VLT, 
i.e. with S < 15°, truncated at z = 0.05. These have been colour coded blue/red for Sy 1/Sy2 where such classifications are available (note 
that only ‘class 4’ AGN are within our luminosity and redshift range). The dotted line is the flux limit for 90% of the sky. The dashed 
lines indicate the distance and luminosity thresholds for our selection in order to have a complete sample - for which the selected targets 
are marked with green. Luminosities in this plot are based on distances derived from the uncorrected redshift, as used in the initial target 
selection. These may differ slightly from the distances listed in Tabled which are, where possible, redshift independent or corrected for 
peculiar velocities. 


focusses on the central 100 pc, assessing whether there 
has been recent star formation, and if so whether it is 
still on-going or has ceased. 

2) What mechanisms are responsible for driving gas 
from the kiloparsec scale into the central tens of parsecs, 
and what triggers these? This includes dynamical pro¬ 
cesses in the central kiloparsec as well as the role of the 
larger scale host galaxy and the influence of the local 
environment. 

3) Do luminous Seyferts always drive ionised and/or 
molecular outflows, and how do they interact with the 
interstellar medium? This addresses the ubiquity and ef¬ 
ficacy of a variety of emission lines as tracers of outflows, 
as well as quantifying the outflow velocities, rates, and 
mass loading. 

These - and other - questions will be addressed 
in future papers using a combination of high spectral 
resolution data fro m 0.4-2.3/im taken with XShooter 
(|Vernet et all 1 20 lH ) and high spatial resolution integral 
field sp ectroscopy in the H- and K-bands taken wit h SIN- 
FONI (|Eisenhauer et al.ll2003t [Bonnet et al .1120041) . both 
on the VLT. At the time of writing, XShooter observa¬ 
tions have been completed for 8 Seyferts, and 10 inactive 
galaxies, and the remainder of the targets are scheduled; 
SINFONI observations are completed or scheduled for 
the first half of the sample. 

2.1. Active Galaxy Sample 

Our first criterion for selecting a sample, and the only 
astrophysical one for the active sample, is to define a set 
of local AGN in a well-characterised way so that the se¬ 
lection effects are (as far as possible) understood. The 
largest optical or infrared surveys that select AGN are 


often incomplete for nearby galaxies, and also tend to 
use AGN tracers that are anisotropic or impose a bias 
against star formation (by using large apertures to mea¬ 
sure features that can be produced by both AGN and star 
formation). For example, [O III] A5007A, a typical optical 
line used in the selection of AGN, is not isotropic: in com¬ 
parison to [OIV]A25.9/rm, which correlates well with 10- 
200 keV emission, it is underluminous in some Sy 2s - es¬ 
pecially those that are Compton thick - b y up to an order 
of magnitude compared to S v Is (IM elende z et alJ 120081 : 
IDiamond-Stanic et alJ 120091 : IWeaver et al.l 1201(1) . The 
same effect is seen in the 2 lOkeV emission when com¬ 
pared to the [OIV] line. However, high ionisation lines 
may not be an ide al way to se lect AGN since they are not 
always observed: iKoss et ahl (]2013D detected NeV only 
in 2/3 of luminous infrared galaxies for which an AGN 
was confirmed by the Swift-BAT survey. 

In contrast to emission line and infrared tracers, the 
very hard 14-195 keV band of the Swift-BAT survey mea¬ 
sures direct emission from the AGN rather than scat¬ 
tered or reprocessed emission, and is much less sensi¬ 
tive to obscuration in the line-of-sight than softer X- 
ray or optical wavelengths (selecting only against highly 
Compton thick AGN). Indeed it is widely accepted as 
the least biased survey for AGN with respect to host 
galaxy properties, and as such it has been well stud¬ 
ied. There is a vast amount of anci llary data o n the 


larger scale host galaxy 

rroDerties (IKoss et al. 

2010 

20111; Winter et al. 2009, 

20101; Melendez et al. 

2014; 

Mushotzkv et al.l 20141). as 

well as analysis of the bolo- 


metric corrections which enable one to m ake an esti¬ 
mate of the AGN luminosities dVasudevan et alJ 120101 : 
IWinter et al.1120 ll : IMelendez et al.ll2014fl . Because it is 


an all-sky survey with roughly uniform sensitivity, we can 
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TABLE 1 

Summary of AGN and host galaxy properties 


Name 

Distance 

AGN 

log h H 

log Li4_i95 

log L 2 _io 

log 

log N h 

Hubble 

axis 

bar 


(Mpc) 

classification 

L 0 

ergs -1 

ergs -1 

ergs -1 

cm -2 

stage 

ratio 


(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 


NGC 1365 

18 

Sy 1.8 a 

10.58 

42.39 

41.83 

42.54 

22 .2 d 

3 

0.55 

B 

MCG-05-14-012 

41 

Sy 1.0 a 

9.60 

42.60 

41.63 


<21.9 

-1 

0.86 


NGC 2110 

34 

Sy 2 (lh) a 

10.44 

43.64 

42.53 

43.04 

23.0 e 

-3 

0.74 

AB 

NGC 2992 

36 

Sy 1.8 a 

10.31 

42.62 

42.05 

42.87 

21.7 

1 

0.30 


MCG-05-23-016 

35 

Sy 1.9 a 

9.94 

43.47 

43.11 

43.42 

22.2 

-1 

0.45 


NGC 3081 

34 

Sy 2 (lh) a 

10.15 

43.06 

41.54 

42.75 

23.9 

0 

0.77 

AB 

NGC 3783 

38 

Sy 1.2 a 

10.29 

43.49 

43.12 

43.47 

20.5 

1.5 

0.89 

B 

NGC 4235 

37 

Sy 1.2 

10.43 

42.72 

41.66 

42.17 

21.3 

1 

0.22 


NGC 4388 

39 

Sy 2 (lh) 

10.65 

43.70 

42.57 

42.93 

23.5 d 

3 

0.18 

B 

NGC 4593 

37 

Sy 1.0-1.2 a 

10.59 

43.16 

42.77 

42.97 

<19.2 

3 

0.74 

B 

NGC 5128 (Cen A) 

3.8 

Sy 2 

10.22 

42.38 

41.50 

41.82 

23.l d 

-2 

0.78 


ESO 021-G004 

39 

Sy b 

10.53 

42.49 

41.21 


23.8 

-0.4 

0.45 


MCG-06-30-015 

27 

Sy 1.2 

9.59 

42.74 

42.51 

42.87 

20.9 

-5 

0.60 


NGC 5506 

27 

Sy 2 (li) 

10.09 

43.32 

42.91 

43.28 

22.4 

1 

0.23 


NGC 5728 

39 

Sy 2 

10.56 

43.21 

41.41 

42.35 

24.2 

1 

0.57 

B 

ESO 137-G034 

35 

Sy 2 

10.44 

42.62 

40.86 


24.3 

0 

0.79 

AB 

NGC 6814 

23 

Sy 1.5 

10.31 

42.69 

42.17 

42.18 

21.0 

4 

0.93 

AB 

NGC 7172 

37 

Sy 2 (li) c 

10.43 

43.45 

42.53 

42.88 

22.9 

1.4 

0.56 


NGC 7213 

25 

Sy 1 

10.62 

42.50 

41.95 

42.58 

<20.4 

1 

0.90 


NGC 7582 

22 

Sy 2 (li)® 

10.38 

42.67 

41.12 

42.81 

24.2 d 

2 

0.42 

B 


Note. — Columns: (1) common name; (2 ) distances are taken from NED, using redshift independent estim ates or peculiar ve locity 
corrections derived by ITheureau et al.| d2007T) where possible (except for NGC5128 where the distance is from [Harris et al.l |201Qt) , and 
so may differ slightly from the distance based on the uncorrected redshift used in the initial selection; (3) AGN type as given in NED, 
except where indicated for the 8 AGN that we can independently classify, with additional information: l i — near-infrared bro ad lines, lh 
= polarised broad line emission; (4) integrated H-band luminosity (given by the 2MASS t otal magnitud ejJSkrutskiee t al.||2006l) as a proxy 
for stellar mass; (5) observed 1 4—195 keV lumino sity (70 month average) from Swift- BAT _|Baum gartner et al.||2013[) ; (6)~observed 2-10 keV 
luminosity (single epoch) from [Ricci et al ] ( 2015); (7) nuclear 12 fim luminosity from [Asmus et^dTif2014l) ('MCG-05-14-012. ESO 021-G004, 
and ESO 137-G034 are not included in this catalogue and, as far as we know, do not have 12 pm measurements on arcsec scales); (8) neutral 
absorbing column, from [Ricci et al ] J2015I) based on modelling 0.3-150 keV spectrum; (9) Hubble stage (NED); (10) axis ratio (NED); (11) 
B and AB denote presence of a bar or weak bar respectively (NED). 

a Confirmed or updated based on our available XSHOOTER data (Schnorr-Miiller et al. in prep.). Note that NGC4593 is ambiguous 
because it is close to the boundary between Sy 1 and Sy 1.2; MCG-05-23-016 has weak broad Ha; NGC 2110 is changed from ‘Sy 2 (li)’ to a 
simple ‘Sy2’, because the polarised broad line has a double peaked profile which is not seen at near-infrared wavelengths, 
k ES0 021-G004 has no optical classification, and our XSHOOTER data for this object have not yet been taken. 

c We have added a ‘li’ qualifier to the classifications of NGC 7172 because lSmajic et al ] <20T2l) show evidence for broad Paa and Br 7 , and 
NGC 7582 because iReunanen et al ] (120030 find broad Br 7 . 

^ Variations in Njy between Compton thick and thin have been reported for NGC 1365 JRisaliti et al.||2009Q and NGC 7582 (Bianc hi et al.| 
I2009Q . Compton thin variations in N# have been reported for NGC 4388 JElvis et al.|[2004Q and NGC 5128 d Beckmann et al.l 1 20 ill) . 
e The absorbing column for NGC 2110 includes partial covering absorbers, and the value given here is weighted by covering factor. Details 
are given in I Ricci et al.| d2015P ; see also lEvans et al.l 120071) . 


select a complete, volume limited sample of AGN. Fur¬ 
thermore, the continuous nature of the survey means that 
more recent versions of t he catalogue (58 and 70 month, 
IBaumgartner et alll2013fl average over variability during 
the last 5-6 years. 

In addition to the single astrophysical criterion de¬ 
scribed above, we also impose two observational criteria 
to ensure homogeneous high resolution and observabil¬ 
ity. Our sample therefore consists of all 20 AGN in the 
58 month Swift-BAT catalogue that meet the following 3 
criteria: (i) 14 195 keV luminosities log Al4—195 > 42.5 
(using redshift distance), (ii) redshift z < 0.01 (corre¬ 
sponding to a distance of < 40Mpc), and (iii) observable 
from the VLT (5 < 15°) so that they tend to be in the 
southern sky. Note that the first two criteria are adjusted 
to intersect the flux limit of the 58-month catalogue for 
90% of the sky, that is 1.46 x 10 -11 ergs -1 cm -2 . The se¬ 
lection, and its completeness, is graphically represented 
in Fig. [T] 


2.2. Matched Inactive Sample 


The second criterion for our selection is that there 
should be a properly matched sample of inactive galax¬ 
ies, so that it is possible to discern which features (stellar 
age, inflow rates, etc.) are related specifically to the AGN 
activity. We describe the selection of the inactive sample 
below. Although it has no part in the analysis presented 
in this paper in Sec. 0] [5J it does play a key role in many 
of the analyses which we will present in future papers. 

The characteristics on which the inactive sample 
matching is based are: host galaxy morphology (Hub¬ 
ble type), inclination (axis ratio) and H-band luminosity 
(as a proxy for stellar mass). For the purposes of target 
selection, we have derived the inclination from the axis 
ratio without any compensation for finite thickness. This 
is not expected to yield a significant bias when comparing 
two galaxies with similar morphological type; and the in¬ 
clinations of the host galaxies will be assessed more care¬ 
fully in future analyses. The relation between H-band 
luminosity and stellar mass has been calibrated from a 
sample of similar galaxies for which the stellar mass has 
been derived using multi-band photometry dKoss et al.l 
IM3). Fig. [2] shows that the scatter in this relation is 
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Fig. 2.— Comparison of the stellar mass M s tar , derived from 
multi-band photometry, to the integrated H-band luminosity Lh 
for a sample of AGN and other galaxies similar to that presented 
here dKoss et al.11201 lT) . There is a clear relation with a scatter of 
0.2 dex, showing that Lh can be used as a proxy for M s tar- 


TABLE 2 

Summary of control sample properties 


Name 

(i) 

Distance 

(Mpc) 

(2) 

log L H 
Lq 
(3) 

Hubble 

stage 

(4) 

axis 

ratio 

(5) 

bar 

(6) 

NGC 718 

23 

9.89 

1 

0.87 

AB 

NGC 1079 

19 

9.91 

0 

0.60 

AB 

NGC 1315 

21 

9.47 

-1 

0.89 

B 

NGC 1947 

19 

10.07 

-3 

0.87 


ESO 208-G021 

17 

10.88 

-3 

0.70 

AB 

NGC 2775 

21 

10.45 

2 

0.77 


NGC 3175 

14 

9.84 

1 

0.26 

AB 

NGC 3351 

11 

10.07 

3 

0.93 

B 

ESO 093-G003 

22 

9.86 

0.3 

0.60 

AB 

NGC 3717 

24 

10.39 

3 

0.18 


NGC 3749 

42 

10.40 

1 

0.25 


NGC 4224 

41 

10.48 

1 

0.35 


NGC 4254 

15 

10.22 

5 

0.87 


NGC 4260 

31 

10.25 

1 

0.31 

B 

NGC 5037 

35 

10.30 

1 

0.32 


NGC 5845 

25 

10.46 

-4.6 

0.63 


NGC 5921 

21 

10.08 

4 

0.82 

B 

IC 4653 

26 

9.48 

-0.5 

0.63 

B 

NGC 7727 

26 

10.41 

1 

0.74 

AB 


Note. — Columns: (1) common name; (2) distances are, where 
possible, redshift independent (NED); (3) integrated H-band lu¬ 
minosity (as a proxy for stellar mass). (4) Hubble stage (NED); 
(5) axis ratio (NED); (6) B and AB denote presence of a bar or 
weak bar respectively (NED). 


0.2 dex. 

We have quantified the matching criteria above using 
a x 2 metric with tolerances of ±1 for Hubble type, ±15° 
for inclination, and ±0.3 dex for H-band luminosity. In 
addition we require that the presence/absence of a bar is 
matched if possible. Finally, also where possible, we have 
selected inactive galaxies with redshifts less than or equal 
to those of their active pairs. All in active g ala xies were 
selected from the RC3 catalogue (Ide Vaucouleurs et al.1 


[199 If) , rejecting any listed in the IVeron-Cettv fe Veronl 
(|2010H catalogue of known AGN. While some may still 
host very low luminosity AGN, the decisive factor is that 
(based on the lack of any obvious signature of BH acre- 
tion in optical, radio, and X-ray data) it is orders of 
magnitude weaker than in the X-ray selected sample. 

The inactive galaxy selection is based on matching the 
characteristics of individual galaxies in the active sample. 
This is used as a robust way to ensure that the samples, 
as a whole, match; it is not specifically our intention 
to compare pairs of active and inactive galaxies. Part 
of the reason is that some inactive galaxies are a good 
match for several active galaxies, which means there is no 
unique pair matching. The final matched inactive sample 
contains 19 galaxies, which are listed in Table [2] together 
with the host properties used for the selection process. 

A graphical comparison of the main host galaxy prop¬ 
erties is given in Fig. [3] By design the distributions of 
these properties should be similar, as is apparent from 
the figure. This is confirmed for Lh , axis ratio, and 
Hubble stage, by Kolmogorov-Smirnov (KS) tests which 
indicate that the differences are not significant. Formally, 
the probabilities that the differences between the active 
and inactive samples could arise by chance exceed 20%. 
Both the active and inactive samples show the same dis¬ 
tribution in H-band luminosity: the active sample has a 
mean of log Lh[Lq] = 10.3 with a ler spread of 0.3 while 
the inactive sample has a mean of 10.2 and a spread 
of 0.4. Both samples cover the full range of axis ratio 
(or equivalently inclination), with about half the sample 
each side of 0.7, equivalent to 45° for a disk. And both 
samples cover a wide range of morphological types, with 
a clear preference for Hubble stage at ~ 1 corresponding 
to SO and Sa ‘early disk’ types. This distributio n incon¬ 
sistent with the analysis presented bv lKoss et al.1 (|2011lh 
who found for a larger sample of Swift-BAT AGN that 
only ~ 10% were in ellipticals, while ~ 30% had inter¬ 
mediate ( SO) and ~ 4 0% ea rly type spiral (Sa-Sb) host 
galaxies. iKoss et all (120111) noted that these fractions 
differ from the distribution of normal galaxies - 16% el¬ 
lipticals, 26% spirals. But the distribution of host types 
of the Swift-BA T AGN are consistent with this once one 
takes into account the detection r ate of Seyferts (rather 
than LINERs ) reported by iHol (120081) , which peaks for 
SO-Sb types. iDavies et al.l ( 20111 ) discussed the role of 
host type and environment in the context of the origin 
of the gas fuelling the AGN, and argued that this has 
an impact on what gas inflow mechanisms one may see 
in the circumnuclear region. This issue of morphology, 
including the presence or absence of a bar, will be revis¬ 
ited for our complete sample in a future work. The main 
point here is that the difference in distributions under¬ 
lines the need to include morphology as a criterion for 
matching the inactive sample as we have done. 

The largest difference between the samples occurs for 
the distance distribution which a KS test indicates may 
be marginally significant at a level of 2.5er. This is a 
direct result of our requirement to select, whenever pos¬ 
sible, inactive galaxies that are closer, rather than more 
distant, to their active pair. We note that the active 
and inactive samples have mean distances of 31 Mpc and 
24Mpc respectively, and in both cases the la spread is 
about 9 Mpc. 
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Fig. 3.— Comparison of several host galaxy properties between the hard X-ray selected AGN (orange; solid lines or plusses), and the 
matched inactive sample (dark green; dashed lines or crosses): H-band luminosity as a proxy for stellar mass (far left), axis ratio (centre 
left), Hubble type (centre right), and H-band luminosity vs distance (far right). There is no obvious relation between host galaxy luminosity 
and distance; and as expected by design, the distributions of the inactive sample properties match those of the AGN very well. The biggest 
discrepancy is in distance, and this results from our preference to select closer, rather than more distant, inactive galaxies when possible. 
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Fig. 4.— Comparison of nuclear properties of Sy Is (blue) and Sy2s (red). There are no significant differences between Sy Is and Sy2s 
for any of the luminosity distributions shown. This holds even for the absorption corrected 14-195keV luminosity (not shown), despite 
the tantalising hint of a difference. In contrast, the difference in log Njj is 3.7o\ Far left: the sample was selected according to 14-195 keV 
luminosity, and shows no significant bias towards either Sy Is or Sy2s. Centre left: Given the small numbers of objects, there is also no 
significant difference in the observed 2-10 keV luminosities. Center right: in contrast, and as expected, Sy 2s do exhibit higher absorbing 
columns, with a median of > 10 23 cm -2 in comparison to ~ 10 21 cm -2 for the Syls (note we have marked AGN with no measurable 
neutral absorption as limits). Far right: the nuclear mid-infrared luminosity shows a remarkably narrow distribution, similar to that of the 
14—195 keV luminosity, with no clear difference between Sy Is and Sy 2s. 


3. AGN PROPERTIES 

In this section we compare various properties of the 
Syls and Sy2s in the active sample, as listed in Ta¬ 
ble H] focussing on the absorbing column and its impact. 
We will, as our data become available, update the clas¬ 
sifications of the AGN if necessary. Similarly, and as 
applicable to the specific analyses we will perform, we 
will re-assess the host galaxy properties of the active and 
inactive samples: morphological classifications, whether 
they are barred, and their inclinations. 

The sample was selected purely on £ 14 - 195 . Adopt¬ 
ing a simple c onversion to es timate the bolometric AGN 
luminosity as dWinter et al.ll2012l ) 

log Lagn = 1.12 log £ 14-195 —4.23 

(i.e. bolometric corrections of 5-10 for the luminosity 
range here) and without applying any absorption cor¬ 
rection (see below), our complete volume limited sample 
has a median luminosity of log Lagn = 43.5 with a ler 
distribution of 0.5 dex. 

Table |T| shows that our selection has yielded similar 
numbers of Sy Is and Sy 2s, and the first panel of Fig. [4] 
indicates that any difference in their £ 14-195 distribu¬ 


tions is not significant. This suggests that, for the num¬ 
ber of objects in the sample, any preference to select 
either Syls or Sy2s is not significant. For both £ 14-195 
and a KS test indicates that the probability of 

the difference between the Sy 1 and Sy 2 distributions 
arising by chance exceeds 45%. The difference for the 
observed £ 2-10 is also not enough to be significant. It 
does, however, have a measurable impact on the ratio 
of the observed 14-195 keV and 2 lOkeV luminosities, 
which is 4 for the Sy Is and 17 for the Sy 2s. The same 
characteristic is clearly seen in the left panel of Fig. [5] 
where £ 14-195 is plotted versus observed £ 2 - 10 - Here 
Sy ls/2s are represented by blue/red colour respectively. 
The Sy Is follow the relation between £ 14-195 and L^ 10 
derived by iWinter et al.1 (120091 ) 

log £|To = 1-06 log £ 14-195 - 3.08, 

(in the luminosity range here, this relation is approxi¬ 
mately equivalent to log L‘^ 10 ~ log £ 14-195 — 0.5) while 
many of the Sy2s do not. The reason is due to their 
different absorbing columns Nh , which is denoted by 
open/filled symbols in Fig. [5] And indeed, the observed 
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Fig. 5.— Comparison of nuclear X-ray and mid-infrared luminosities of Sy Is and Sy2s (denoted by blue and red colour respectively). 
Absorbing column is represented by open/filled symbols: filled squares have Nh > 10 23 cm -2 while open squares have Nh < 10 23 cm -2 
(as discussed in Sec. \3\ this thresh old is where Nh starts to have an impact on the 2-10keV band). The two red open squares are both 
AGN classified as Sy2(li). IBurtscher et all (1201511 argue that, based on their physical properties, these should be considered more like 
Sy 1 than Sy2. Left: Li 4 _i 95 vs observed L 2 — io- From top left to bottom right, the dotted lines trace ratios of L 14 — 195 /L 2—10 = 1, 10? 
100 (indicative of increasing absorbing column, see lKoss et al.lf2013l) . The relation of I Winter et al.l (|20090 for is indicated by the 

dashed line. Sy Is and other AGN with low N# follow this relation, indicating that for these L 2—10 ~ ^ 2 ^ 10 ■ Right: ^ 14—195 vs nuclear 
Li 2 fim• From top left to bottom right, the dotted lines indicate ratios of Z/ 14 —i 95 /Z/i 2 £ 4 m — 1? 10- The dashed line denotes the relation of 
I Gandhi et al.l <120100 , replacing 10 with L 14—195 according to I Winter et alj (I2009Q . To first order, the data follow this relation without 
any dependency on Nh or AGN type. 


^14-195/^2-10 r atio can in princ iple be used to give an 
estimate of Nh ( Koss et aPf2013 i. 

The difference in Nh is clearly seen in the centre right 
panel of Fig. []] and a KS test gives it a significance o f 
3.7cr. Values for Nh are taken from iRicci et al.l (|2015H . 
who derived absorbing columns for many of the Swift- 
BAT AGN in a consistent way, fitting the 0.3-150 keV 
spectral energy distribution. The basic model includes 
an absorbed cutoff power-law continuum plujy a reflec¬ 
tion component (for details see IRicci et al.l [20T5h . For 
sources with Nh < 10 22 cm -2 , they included various ad¬ 
ditional components if statistically needed: a blackbody 
(for the soft excess), partially covering ionized absorp¬ 
tion, a cross-calibration constant (for possible variabil¬ 
ity between the non-simultaneous soft X-ray and Swift- 
BAT observations), and an iron line (or other emission 
lines in that region). Similarly, for the more obscured 
sources, some additional components were added if re¬ 
quired: a scattered component, a collisional plasma, 
emission lines, and a cross-calibration constant. Note 
that in the column densities given in Table[l]we do not in¬ 
clude warm absorbers, which arise in ionised o ut-flo wing 
gas and are observed in about 50% of Sy ls dKomossal 
119991 : iBlustin et, al.ll2005l : iWinter et, al .1120121) . The differ¬ 
ence between the neutral absorption in Sy Is and 2s is sig¬ 
nificant. The mean for the Syls is log-/V#[cm -2 ] ~ 21.1 
while for Sy2s it is log Nh [cm~ 2 ] = 23.5. 

Models of the propagation of X-rays through 
gas in various geometries around a centra l source 
dMatt et al.l Il999t iBrightman fe Nandral 1201 It lYaaoobl 
20121 IRicci et al .1120151) indicate that Nh can account for 
the different luminosity ratio between the Sy Is and Sy 2s. 
These models indicate that the flux correction for the 14- 
195 keV band is negligible up to 10 23 cm -2 , and reaches 
a factor 2 around 2-3 x 10 24 cm -2 , and a factor 10 at 


10 25 cm' 2 (whether higher corrections for larger columns 
are required depends on the reflection component). As 
such, even the absorption corrected 14-195 keV luminosi¬ 
ties do not, based on a KS test, yield a significant differ¬ 
ence between Sy Is and Sy 2s. On the other hand, large 
flux corrections, greater than a factor 2, to T 2-10 are 
already required for absorbing columns of ~ 10 23 cm -2 , 
increasing to a factor 10 at ~ 5 x 10 23 cm” 2 . If not 
accounted for, this can lead to a bias when selecting via 
the 2-10 keV band, since the impact on the observed flux 
means that some Sy 2s may be excluded from a sample 
even before such a correction can be made. 

This effect is much reduced when selecting in the 14- 
195 keV band. Its impact is apparent in t he distribu- 
tion of N h for the Sy 2s: as expected from iRisaliti et ahl 
(HU, most of the sources have N h > 10 23 cm 2 . And 
(allowing for variations in N h) 3-5 of the sources have 
N h > 10 24 cm^ 2 , consistent with t he expectation that 
10-30% of AGN are C ompton thick (iTreister et al.ll2009l: 
lAlexander et al.l I2013T) . But in contrast to IRisaliti et al l 
(j!999f) there are none with N h > 10 25 cm -2 , which is 
simply due to the absorption in the band at such high 
columns. However, we expect very few local sources to 
fall in this regime. NGC 1068 is the best known example 
of one that does, and hence is excluded from the sam¬ 
ple: although its ^ 14 - 195 /^ 2-10 ~ 10 is not im mediately 
indic ative of a high column, early estimates (IMatt et al.l 
119971) suggested it is highly Compton thick, and fits to 
recent hard X-ray data requi re sever al components with 
columns up to 10 25 cm -2 dBauer et al.ll2014l ). 

The final compari son is o f the nuclear 12/rm lumi¬ 
nosities, taken from lAsmus et all (j2014f ). These data 
were observed at subarcsecond resolution, and so re¬ 
solve out mid-infrared continuu m due to sta r formation 
in the circumnuclear region. iGandhi et al.1 (|2010T) and 
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lAsmus et al.l (l2012f) have shown that there is a very good 
correlation between L l ff 10 an d with no major dif¬ 

ference between Sy Is and Sy 2s. Our comparison in 
Figs. [4] and [5] confirms that this is also the case for our 
sample, that there is little dependence on or its 

relation t o Am-iar with e i ther A GN type or Nh- Note, 
however, iBurtscher et al.l (12015T) point out that the nu¬ 
clear 12/iin continuum does show some signs of a slight 
anisotropy, which can also be seen in the right panel of 
Fig. [5] 

4. OPTICAL OBSCURATION AT LOW LUMINOSITY 

The proportion of AGN with substantial optical ob¬ 
scuration and/or X-ray absorption, and their depen¬ 
dence on luminosity, offers insight into the properties of 
the obscuring structures around the AGN. In this sec¬ 
tion we focus on low luminosities, in the range 42.5 < 
log -Li 4 _i 95 < 44. After clarifying in Sec. 14.11 the defini¬ 
tions we use for optically obscured and X-ray absorbed 
AGN, and summarising the various types of AGN that 
are found, we discuss X-ray unabsorbed Sy 2s in Sec. 14.21 
and put this in the context of the Swift-BAT sample in 
Sec. 14.31 to assess how rare such objects really are. 

4.1. Definitions and Classifications 

Much of the discu ssion below is set in the context of 
iMerloni et al.1 (I2014T) . As such, we follow these authors 
and adopt N h > 10 215 cm -2 as the criterion for defining 
an AGN to be X-ray absorbed. For typical gas-to-dust 
ratios this is equivalent to an optical screen extinction 
of A y ~ 2 mag, sufficient to mildy obscure the line of 
sight to the broad line region (BLR) at visible wave¬ 
length^. In our volume limited 14-195 keV sample, we 
find 50-60% of the AGN are optically obscured while 
~ 60% are X-ray absorbed. Most of the obscured AGN 
are both optically obscured and X-ray absorbed; none 
are just optically obscured; and only three are just X- 
ray absorbed. These exceptions are ‘minor’ since they 
are all close to the boundary between the regimes. The 
X-ray absorbed AGN NGC1365 is optically classified as 
Syl.8 and, although in many studies is taken to be a 
Sy2, here we adopt a strict definition and so classify 
it as Syl. Similarly MCG-05-23-016 is X-ray absorbed 
and has only weak broad Ha, hence its previous classi¬ 
fication as a Sy2 which we have revised to Syl.9. The 
Sy 1.8 NGC 2992 has a column only marginally above our 
X-ray absorption threshold and so is on the borderline 
in both classifications. It would be classified as X-ray 
unabsorbed according to other definitions that se t the 
threshold at 10 22 cm -2 ( Panessa fe Bassanl 120021) . As 
such, there could be two more optically obscured and 
one less X-ray absorbed AGN than found using the strict 
definitions we apply in Fig. [6] This is reflected by the 
ranges shown in Fig. [7] 

The joint optical/X-ray classifications are sum¬ 
marised in Fig. [6] where we adopt the terminology of 

b We do not distinguish a category of obscuration in which the 
broad emission lines are obscured at optical wavelengths but can 
be detected at near-infrared wavelengths. These objects, which 
are sometimes classified as Syli, are considered to be Sy2s in 
this work. We note, however, that like Syl.8-1.9s they are close 
to the boundary of optically obscured/unobscured and X-ray ab¬ 
sorbed/unabsorbed. And there are good reasons to co nsider them 
instead as Sy 1 as discussed by IBurtscher et al] 1120151 1. 


42.5 < log L 14 _ 1g5 [erg s '] < 43.6 



Fig. 6. — Fractions of optically obscured and X-ray absorbed 
AGN in our complete volume limited Swift-BAT sample (see also 
Fig-0- The first/second digit of the type codes for optical/X-ray 
obscuration as described in Sec. 14.11 The black lines denote the 
strict definitions. The grey lines show the impact of relaxing the 
definitions slightly: the three type 12 AGN would be reclassified 
as one type 11 and two type 22. 



Fig. 7. — Fractions of optically and X-ray obscured AGN as a 
function of lu minosi ty. The curves are adapted from Fig. 12 of 
IMerloni et al.l I12014J) . The details o f th e types, defined by those 
authors, are summarised in Section 14.31 the first/second digit of 
the type codes for optical/X-ray obscuration. The luminosity scale 
has been derived from as indicated in Section [3] The black 

ranges refer to our sample, and indicate the impact of allowing 
some flexibility in the definition of optically obscured and X-ray 
absorbed. At overlapping luminosities, the sample have similar 
fractions of AGN that are X-ray absorbed (types 22+12) and that 
are both X-ray absorbed and optically obscured (type 22); but the 
fractions of AGN classifi ed as Sy 2 (types 22+21) are very different. 
IMerloni et al.l (l2014ll in which the first/second digit of 
the classification denotes the optical/X-ray type. This 
leads to the following types: 

type 11: optically unobscured and X-ray unabsorbed; 

type 22: optically classified as Sy 2 and X-ray absorbed; 

type 21: optically classified as Sy2 but X-ray unab¬ 
sorbed; 

type 12: optically unobscured but X-ray absorbed; 

Thus, with ref erence to the curves in Fig [7] adapted from 
IMerloni et al.l (I2014T) to the 14-195 keV luminosity scale, 
we also have: 
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type 22+12: all X-ray absorbed AGN (regardless of 
their optical classification); 

type 22+21: all AGN optically classified as Sy2 (inde¬ 
pendent of whether they are X-ray absorbed). 


Comparing the AGN in our sample to the left edge 
of Fig [7] immediately highlights a discrepancy. Above, 
we noted that 50-60% of our complete sample are opti¬ 
cally obscured and cl assed as Sy2s, a f r action similar 
to that esti ma ted bv [Lawrence fe Elvisl (120101 ). How¬ 
ever, iMerloni et al.l ([20141 1 class ~ 90% of their AGN 
in an overlapping luminosity range log £ 14-195 ~ 43- 
43.5 as Sy2s. The difference is due to the large frac¬ 
tion of AGN classified by these autho rs as type 21 , that 
is X-ray unabso rbed _Sy 2 galaxies dPappa et all 
Panessa i- & - ^assajii l2002t IBrightman fe Nandral 


Bianchi et al.112012 1. 


2001 


2008 


4.2. How common are X-ray unabsorbed Seyfert 2s? 

X-ray unabsorbed Sy2s are believed to be galaxies 
with a direct view to the AGN but in which there is 
no BLR (such objects are also known as pure or true 
Sy2s). This view is supported by, for example, the 
6 objects identified by lHawkinsI (l2004lj as having opti¬ 
cal spectra typical of Sy2s (H/3 FWHM < lOOOkms -1 
and [O III]/H/3 > 3) but large amplitude variations typ¬ 
ical of Sy Is. Subse quent X-ray observations of three by 
IGliozzi et al.l (1200711 confirmed that these do not have sig¬ 
nificant absorption in the 0.3-8 keV band. These authors 
also showed that these AGN are relatively luminous (in¬ 
trinsic log £ 0 . 5-8 > 43.2). Thus, despite t heir rath er low 
Eddington ratios Lboi/LEdd < 0.01 ( Gliozzi et al.ll2007lf . 
the absence of a BLR cannot be explained by low lu¬ 
minosities and/or ac cretion rates ei t her being unable to 
sustain a disk wind ([Nicastrol 1200(1 : lElitzur fc Hoi 12009]) 
or meaning that high dispersion prevents BLR clouds 
from surviving (lLaorll2003h . At present, the reason that 
they may not have a BLR is still open. 

These objects may be related to the Sy2s for which 
spectropolarimetry shows no evidence of a hidden BLR. 
There have been many studies of polarised emission from 
hidden BLRs. The largest indicate that Sy2s in which 
a hidden BLR has not been detected have less luminous 
AGN t han other Sy2s but are not more optic ally ob¬ 
scured (|Tranll20f)illGu fc Huandl2002tlTranll2003H - qual¬ 
itatively, but not necessarily quantitatively, consistent 
with the idea that BLR may not form in low luminosity 
AGN. It also suggests that a BLR is completely absent 
in at least 50% of Sy2s. However, this fraction is diffi¬ 
cult to put in perspective because of complications due 
to the impact of detection limits as well as the expected 
scatter in g effici ency. The latter issue was addressed by 
iHeisler. Lumsden. fc Bailevl ( 1997 1 who showed that the 
ability to detect polarised emission was related to the far- 
infrared colours in a way that suggested the scattering 
particles were in the throat of the torus, and so even they 
would b e hidde n at high inclinations. Following on from 
work of lRamos Almeida et all ( 2011 ) which shows that 
the intrinsic properties (rather than just the inclination) 
of the tori in Sy Is and Sy2s are different, Ichikawa et al. 
(120151) have fitted torus models to Sy2s with and with¬ 
out polarised emission. Their results again indicate that 
whether a hidden BLR is seen depends on the location 


of the scattering material in the throat of the torus; but 
that rather than inclination, it is the opening angle and 
covering factor of the torus that determines this differ¬ 
ence. As such, the fraction of Sy 2s in which polarised 
scattered BLR emission is not observed could be very 
different from the fraction without a hidden BLR. Since, 
based on these models, the polarised BLR emission may 
in some cases be only relatively lightly obscured, deep 
and/or near-infrared spectropolarimetry may shed fur¬ 
ther light on this issue. NGC 7 172 and NGC 7582 (see 
Tabled]) are cas e s in p oint here: iBian fc Gul (|2007t ) and 
iMarinucci et alJ ([ 2012 ) both repo rt that no po l arised 
broad Ha has been detected but iSmaiic et alJ (| 2012 H 
have seen broad Paa a nd Bry directly in NGC 7172 and 
iReunanen et al.l (120031) broad Bry in NGC 7582. 

Despite the uncertainties, the evidence that at least 
some relatively luminous AGN do not have a BLR is 
convincing. The key question for this paper is how com¬ 
mon such objects are, since this can have a significant 
impact on the total fraction of AGN optically classified 
as Sy2. 

In their analysis. iRisaliti et al.l (1 9991 ) found that about 
4% of Sy2 galaxie s had Nh < 10 22 cm -2 . However, 
iPanessa fe Bassanil (120021) presented a sample of 17 ob¬ 
jects which had been classified as Sy2 in the literature, 
and which had similarly low column densities. They sug¬ 
gested that X-ray unabsorbed Sy2s may be much more 
common, being 10-30% of all Sy2s. More recently, other 
authors have confirmed that such objects do exist, but 
without a consensus on how common they are, because 
many candida tes are rul ed out on cl oser examination. 
For example, IBrightman fc Nandral (I2008D found that 
of 8 candidates, 4 were underluminous in the 2-10 keV 
bands suggesting that they were Compton thick - with 
the unabsorbed softer X-rays originating in a scattered 
continuum or from host galaxy contamination. Similarly, 
in a n examination of 8 c a ndida tes - includ ing some from 
t he |Panessa fc Bassanil ( 2002 ) sample - iBianchi et alJ 
(|2012D ruled out 4 for a variety of reasons, and confirmed 
only 3 as unabsorbed Sy 2s. Using a variety of metrics 
for a critical assessment of 24 objects that were previ¬ 
ously claim ed to be unabsorbed Sy2 s (again in c luding 
many from IPanessa fc Bassanil 2002 ). iShi et al.l ( 201 0) 
could confirm only two as genuine, with broad emission 
lines 2-3 orders of magnitude fainter than typical Sy Is - 
although several more candidates may have anomalously 
weak broad lines, this could not be confirmed with ex¬ 
isting data. These authors concluded that 1% or less of 
Sy2s are X-ray unabsorbed. 

In contrast, a much higher fraction of such sources 
- 30% of AGN at l og £14-195 ~ 43 - appears in the 
IMerloni et al.l (|2014D sample. Whether this indicates 
that there is a large population of X-ray unabsorbed Sy 2s 
or is due to a classification bias was discussed in depth by 
these authors. With reference to stacked optical spectra 
and spectral energy distributions, they suggested that 
many of the lower luminosity AGN may have been incor¬ 
rectly photometrically classified as Sy2 due to the low 
contrast of the AGN against the relatively bright host 
galaxy. On the other hand, spectroscopic data are sensi¬ 
tive to features associated with Sy Is, such as broad lines, 
even when they are weak. Such signatures would not 
be evident in broad-band integrated photometry. This 
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difference is reflected in their data: while 77% of the 
AGN with only integrated broad-band photometric data 
are classified as Sy2, only 56% of AGN with spectro¬ 
scopic data are classified as Sy 2s. The lower total frac¬ 
tion of Sy 2s based on spectroscopic classifications implies 
in turn a much lower fraction of X-ray unabsorbed Sy 2s. 

Fig. [ 6 ] shows that in our volume limited sample of Swift- 
13 AT AGN, there are no sources classified as X-ray un¬ 
absorbed Sy2, implying these sources are rare. In this 
context, a simulation of the AGN population can provide 
valuable insight, since we know the Sy2 fractions in both 
flux limited and volume limited samples. We describe 
this simulation in Sec. £3 aimin g to answer the ques¬ 
tion of which curve from IMerloni et all (12014) in Fig. [T] 
properly traces the fraction of Sy 2s: is it types 22+21 or 
just type 22 ? 

4.3. A Simulated Swift-BAT sample 

Our aim here is to simulate parent populations of 
AGN with different prescriptions for the intrinsic frac¬ 
tion of Sy2s. After applying observational limits, we 
can then use c onstr aints from the flux-limited sample of 
IWinter et al.1 (120091 ) and our own volume-limited sam¬ 
ple to discriminate between the different prescriptions, 
giving us a handle on the true fraction of Sy 2s. 

We begin by constructing a population of AGN using 
a 2-10keV luminosity function from [Aird et al.l ( 201 0 s ) 
suitable for low redshift: 

<£(4-io) [(^lo/^o) 070 + (i&o/^o) 314 ]- 1 

where log£o = 44.96. We have then randomly assigned 
each AGN to be a Sy 1 or Sy 2 according to the obscured 
fraction as a function of luminosity. We use 2 functions, 
which are discussed below. We allow the AGN to be 
observed by calculating £44-495 from L l ff 10 , assigning 
them random locations in a local volume, deriving their 
fluxes, and assessing whether they would be detected in 
the Swift-BAT survey. In doing this, we have not in¬ 
cluded the effects of X-ray absorption in the 14-195 keV 
band since, as we have seen, it has little impact at such 
high energies for the absorbing columns expected. We 
use a limit appropriate for the 9-month survey in or¬ 
der to_^elect_aflux limited sample matching that used 
bv IWinter et al.l (120091k and a limit appropriate for the 
58-month survey to match our volume limited selection 
outlined in this paper. 

The left panels in Fig. [ 8 ] show the 2 obscuration func¬ 
tions used in our simulations. The d otted line in the up¬ 
per left panel is taken directly from IMerloni et ahl (|2014f ) 
and was designed to follow the curve for types 22 + 21 ; the 
dashed line in the lower left panel is our modification of 
this function to match the type 22 data, those AGN that 
are both optically obscured and X-ray absorbed. The 
results of our simulations using these functions are also 
given in Fig. [5] The centre panels show the Sy 1 and Sy 2 
populations generated - which can be compared to the 
real data in Fig. [T] The right panels show their resulting 
distributions as a function of luminosity, which can be 
compared to Fig. [9] In these panels, the sharp decrease 
in the number of Sy2s around log £ 44-495 ~ 44 sim¬ 
ply reflects the rapid change in the fraction of obscured 
AGN at this luminosity. The flux limit means that low 
luminosity AGN are not detected at larger distance, and 


hence biasses the apparent fraction of obscured AGN at 
log £44—495 < 43.5. In contrast, because our sample is 
complete to ^£ 44-495 = 42.5, it provides an unbiassed 
indication of the obscured fraction at these luminosities. 

It is immediately clear that the lower panels better 
match the o bserved distr ibution of Sy 1 and Sy 2 sources 
reported in IWinter et al.l (l2009f) . Similarly, the lower 
panels correspond to finding a fraction 54% of Sy 2s (in 
contrast to 87% for the upper panels) in a complete vol¬ 
ume limited sample equivalent to that described in Sec. [3 
again providing a good match to the fraction actually 
measured in the Swift-BAT sample. 

Our conclusion is that, over the luminosity range 43 < 
log £44—495 < 44.5 where the ‘type 22’ and ‘types 22+21’ 
curves differ in Fig. [3 the true fraction of Sy 2s is given 
by the ‘type 22’ curve. This confirms that X-ray un¬ 
absorbed Sy 2s are rare, at most a few percent of the 
Sy2 population; a nd that the high number of type 21 
objects in the IMerloni et al.l (l2014f) sample is due to a 
bias in photometric classification. It also implies that, in 
low luminosity systems, optical obscuration and X-ray 
absorption are usually found together - and hence most 
likely originate in the same shared obscuring structure. 

5 . THE TORUS AND BROAD LINE REGION 

Fig.0shows that the behaviour of the Sy 2 fraction and 
X-ray absorbed fraction of AGN are different above and 
below a luminosity of log £ 44-195 ~ 44. In this section we 
explore what this may indicate about the torus. Often, 
the torus is simply considered to be a dusty entity, since 
it must obscure the BLR in some AGN. A more careful 
definition would be that the torus is a geometrically thick 
structure (since it confines the ionisation cone and any 
outflow) that causes optical obscuration and some (per¬ 
haps much) of the cold X-ray absorption. We show that, 
in this case, its inner boundary is decoupled from the 
dust sublimation r adius. Thus, while near -infrared rever¬ 
beration ma pping rtSuganum a et all 1 20061) and interfero¬ 
metr ic data (jKishimoto et a.1.112009L l20flt iWeigelt et all 
120121 ) measure the dust sublimation radius, they do not 
measure the inner boundary of the torus. 

5 . 1 . At low luminosity, log £14-195 < 44 

At luminosities below log £ 44-195 ~ 44, the optically 
obscured AGN fraction of ~ 60% is similar to the X- 
ray absorbed fraction, and is independent of luminos¬ 
ity. Since the dust sublimation radius Rdust increases 
with luminosity, this independence must reflect the geo¬ 
metrical properties of the torus which we can constrain 
as follows. A luminosity of log £44-195 = 44 implies 
Rdust. 44 ~ 0.2 pc for graphite grains of size 0.05 /im 
(jNetzen I2014D . An obscured fraction of 60% implies 
H/R ~ 0.75 - a half opening angle of ~ 53° - and 
hence Hdust,a ~ 0.15 pc. And that the obscured frac¬ 
tion is constant with luminosity implies that H/R can¬ 
not increase at radii smaller than Rdust, 44 (and the rapid 
drop in the optically obscured fraction at higher lumi¬ 
nosities indicates that H/R also decreases at larger radii 
out to a few times Rdust, 44 )- This is illustrated in the 
left side of Fig. [lOl There are two implications. First, 
it is the same geometrically thick structure that causes 
the optical obscuration and the X-ray absorption, and 
this structure is synonymous with the dusty torus. The 
implied half-opening angle of ~ 53° is at the upper 
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Fig. 8. — Simulations of the Swift-BAT AGN sample using the 2 different optical obscuration fractions, as described in Section 14.31 Top 
panels: using the fraction of AGN classified as Sy 2 (type ‘22+21’); bottom panels: using the true fraction of optically obscured AGN (type 
‘22’). Left panels are from Fig. [7] showing the fraction of AGN classified by I Merloni et ahl (120141) as Sy2 (types 22+21, top) and the Sy2s 
that are also X- ray absor bed (type 22, bottom). The dotted line (top) indicates a functional approximation to the ‘types 22+21’ curve 
from I Merloni et all (|2014l ): the dashed line (bottom) is our modification to match the ‘type 22’ curve. Centre panels show each AGN in the 
luminosity-redshift plane (for clarity only a random subset is drawn). The flux limit is clearly apparent, and the thresholds for our complete 
sample are marked. These panels can be compared to Fig. fTl Right panels: histogram s showing the distrib utions of AGN luminosity in a 
flux limited sample, that can be compared with Fig. [9] which is redrawn from Fig. 7 in I Winter et al.l (120091) . In all panels, Sy Is are drawn 
in blue and Sy 2s in red. 



lo 9 L, 4 _, 95 ke v [erg s _l ] 

Fig. 9.— Numb er of Sy Is (blue) and Sy2s (red) in the flux 
limited sample of lWinter et all l2009i ). This figure is redrawn from 
their Fig. 7 using data from their Tables 1 and 4, but plotted as a 
function of 14 195 keV luminosity derived from ^ 2 -\ 0 as indicated 
in Section [3] The luminosity range to the right of the dotted line 
is that discussed in this paper. 

end of the 30-50° range deri ved from 3D models of the 
ionisation cone kin ematics (Muller-San chez et al.1 120111 : 
iFisclier et al.ll2013f l. a marginal consistency that could 
reflect uncertainties in the derivations of these numbers 
or may have a more physical origin. Secondly, the in¬ 
ner wall of the dusty torus - which defines the optically 
obscured and X-ray absorbed fractions - has universal 
properties at logLi 4 _ig 5 < 44: the average torus prop¬ 
erties are an inner radius of ~ 0.2 pc and a scale height 
of ~ 0.15 pc, independent of luminosity. Here we make 
two notes. The first is that this is a constraint on the 
geometrically thick nature of the torus, and any dust 
that exists at smaller radii will not contribute to the ob¬ 
scured fraction because it is not driving the maximum 


value of H/R dNetzer &; Laod[T~993t iNetzerl 120141) . The 
second is that at luminosities below that of our sample, 
one would expect the torus properties to change if there 
is a lumi nosity t h resh old below which the torus disap¬ 
pears (lElitzur fe Shlosmanll2006l) . 

5.2. At high luminosity, logLi 4 _ 195 > 44 

While none of the AGN in ou r sam ple are this lumi¬ 
nous, the data of lMerloni et al.1 (I2014D . as reproduced in 
Fig. H show a dramatic difference in the Sy2 and X-ray 
absorbed fractions. At these luminosities the Sy2 frac¬ 
tion is not expected to suffer a significant bias due to 
mis-classification simply because the AGN are brighter 
with respect to the host galaxy ; and in deed the fraction 
of AGN classified as type 21 bv lMerloni et al.l (120141) does 
decrease dramatically. As such, one can consider the Sy 2 
fraction as the optically obscured fraction. This fraction 
drops rapidly to ~ 20 % following the L °' 5 dependence of 
the dust sublimatio n radius, as understood in the ‘reced¬ 
ing torus’ concep t (lLawrencelll99H lSimpsonlll998l 120051 : 
iLusso et al.l 12014)1 . without requiring any change in the 
scale height of the torus. In contrast, the X-ray absorbed 
fraction remains nearly constant, falling only from 60% 
to 45%. 

The implicit pop ulation of X-ra y absorbed Sy Is were 
discussed by [Merloni et al.l (I2014D who argued that the 
X-ray absorption is not occuring on large galaxy-wide 
scales and suggested instead that it is more plausibly due 
to dust-free gas within the BLR. Below, we argue that 
this is indeed likely, since BLR clouds are expected to 
contain a significant column of neutral gas; and associate 
the absorption with a ‘neutral torus’ that co-exists with 
the BLR, and which has an impact on how the ‘receding 
torus’ concept should be understood. 
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Fig. 10. — Cartoon illustrating how the inner edge of the torus and the dust sublimation radius Rdust are related, depending on luminosity 
(we do not distinguish smooth or clumpy torus concepts here, and the cartoon could also be envisaged as a less precisely demarcated cloud 
distribution). At the luminosity range on the left side, the fraction of optically obscured and X-ray absorbed AGN is the same. Because 
Rdust is smaller than the radius of the inner edge of the geometrically thick dusty torus, it has no impact on the obscured fraction (which 
is set by the height of the thick structure). The right side depicts the situation at higher luminosities. Here, one finds essentially the same 
geometrically thick gas structure, but now Rdust is larger than the radius of its inner wall. Thus the broad line region extends outward 
into the neutral torus which, further out still, becomes the dusty torus. At these luminosities, the fraction of X-ray absorbed AGN is the 
same, but the optically obscured fraction decreases with increasing luminosity. 


However, first it is important to clarify whether warm 
absor bers, which hav e been observed in about 50% of 
Sy Is (IKomossal 119991: Blustin et al.1 l2005t iWinter et al.1 
Ml ILaha et al .1 12m4fr may be responsible for the X- 
ray absorbed Sy Is here. Warm absorbers imprint their 
signature on an X-ray spectrum through absorption lines 
and are most prominent in soft X-ray bands. They are 
classically identified by the O VII and O VIII absorption 
edges at 0.74 keV and 0.87 keV respectively, ra ther t han 
through their impact on the spectral sha pe (IKomossal 
Il999t [Winter et all 120121 ILaha et all I2014T) . Warm ab¬ 
sorbers have a smaller impact than neutral gas on the 
global sh ape of th e spec trum, as evident from the mod¬ 
els of iPage et al.1 (2 011 ) whose fits to 0.2- lOkeV QSO 
spectra show that a warm absorber requires an order of 
magnitude more gas than a neutral absorber to repro¬ 
duce the same s pect ral shape. The median ionised col¬ 
umn derived bv IWinte r et al.| (j2012) for the Swift-BAT, 
AGN as well as bv ILaha et al.1 ( 2014T) on a different sam¬ 


ple, is ~ 10 2i cm - "; with columns an order of magnitude 
lower for AGN without strong OVII and OVIII detec¬ 
tions. While significant, this is less than ~ 10 215 cnr 2 
threshold adopted here for defining X-ray absorb e d sys - 
tems. Coupled with the fact that iMerloni et ak !j (l2014f ) 
used either the full 0.5-10 keV band or the 0.5-2 keV ver¬ 
sus 2-10 keV hardness ratio to derive the absorbing col¬ 
umn from its effect across the whole band rather than 
in specific features, this strongly argues that the absorp¬ 
tions they derived are due to neutral rather than ionised 
gas. 

The BLR provides an ample source of neutral gas 
that could be responsible fo r the d erived absorption. 
Photoionisation models by iNetzefi (|2013f ) show that 
for a typical BLR ionisat i on parameter of U ~ 0.01 
dLeighlv fc Casebeeil l2007t iNegrete et al.1 1 201. ‘Il l and a 
density of 10 10 cm -3 the temperature drops below 10 4 K 
at a column of ~ 10 21 cm -2 and the fractional abundance 
of HI increases dramatically. Thus, for clouds with a col¬ 
umn density of 10 23 cm” 2 as typically expected, the ma¬ 
jority of gas is neutral rather than ionised. This has been 


used to explain the cause of Sy Is which change state on 
short timescales, in terms of a BLR cloud passing across 
the li ne-o f-sight to the central engine (se e iRisaliti et al.1 
120101 and iTorricelli-Ciamponi et al.ll20l4 and references 
therein). In this interpretation, the clouds have cores 
with column densities of at least a few x 10 23 cm” 2 that 
cause the neutral absorption, and i onised ta il s th at are 
the o rigin of the warm absorption (Risal iti et al.l 12009. 

[Ml . 

Since the BLR clouds contain abundant neutral gas, 
it is natural to associate them with the X-ray absorbed 
Sy Is. The remarkable constancy of the fraction of X- 
ray absorbed AGN across the luminosity range in Fig. [3 
suggests that the geometry of the cold gas absorber also 
changes rather little: both the inner edge and scale height 
of the thick gas distribution remain similar to those at 
low luminosities. The implication is that there is a ‘neu¬ 
tral torus’ that co-exists with the BLR, and extends out 
to the start of the dusty molecular torus as illustrated in 
the r ighthand side of Fig. flOl In this co ntext, w e note 
that iMinezaki fc Matsushital (|2015fl and iGandhi et al.1 
(|2015fl have independently suggested the core of the fluo¬ 
rescent Fe Ka line may originate at radii extending from 
the BLR out to the dusty torus, implying the presence 
of gas at intermediate radii. There is no clear boundary 
between the BLR and the neutral torus since it is the 
same clouds that contribute to both. The BLR emission 
will follow a radial dependence associated with the ioni¬ 
sation parameter. Since U oc Lagn /(-R 2 nu), the rapid 
decrease of U with radius R may be largely mitigated if 
the cloud density nn also decreases with R. Thus the 
BLR could extend far into the neutral torus. How far 
it extends is give n by the scaling o f the BLR size with 
AGN luminosity (Kaspi et al. 20001) . Using the most re - 
cent measurement of that relation dBentz et, al. |[2013). 
and a dopting a ratio Lagn/ AL 5100 ~ 15 (iGrupe et all 
Ml . we estimate the characteristic radius of the BLR 
to be a factor of a few less than the dust sublimation ra¬ 
dius. However, it is also possible that the BLR emission 
could still occur as far as the dust sublimation radius. 
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In the context above, the idea of the ‘receding torus’, 
in which the location of the inner wall of the torus is 
set by the dust sublimation radius, can be misleading. 
We argue that the inner dust boundary and the inner 
gas boundary of the torus should be considered sepa¬ 
rately. The discussion above suggests that the location of 
the inner gas boundary of the geometrically thick struc¬ 
ture (causing the X-ray absorption), is roughly indepen¬ 
dent of AGN luminosity in the range 43 < log £ 14-195 < 
45. And it is within this geometrically thick gas struc¬ 
ture that the inner dust boundary (i.e. the ‘reced¬ 
ing torus’) simply represents the location at which the 
predominant gas phase changes from ionised/neutral to 
dusty/molecular. The thick gas structure itself remains 
nearly unchanged. This means that at low luminosities, 
the thick gas structure is synonymous with the standard 
dusty torus; but at high luminosities, the thick gas struc¬ 
ture has an outer part that is the dusty torus, as well as 
an inner part that contributes to the BLR emission but 
also acts as a n eutral dust-fre e torus. 

Intriguingly, iWinter et all (l2012h found that the 
strength of the O VII and O VIII edges tracing the warm 
absorbers in Swift-BAT AGN are correlated with the 
neutral column density - suggesting that material in the 
neutral torus is associated with a more-or-less propor¬ 
tional amount of ionised material. This can be put in 
the context of AGN showing rapid variations in absorb¬ 
ing column, i n which neutral clouds ha ve an ionised out¬ 
flowing tail (iRisaliti et al.l 12009 . [20111 ) - an interpreta¬ 
tion that applies to individual clouds in the BLR. Given 
that warm absor bers are outflowing dBlustin et al.l120051 : 
lLaha et ~aTTI 1 20 1 111 , and that the distance of the warm ab¬ 
sorbers from the AGN is typically larger than the BLR 
but within a factor of a few of the dust sublimation ra¬ 
dius (|Blustin et al.l 120051) . one might sp e culate that the 
ionised material seen by IWinter et al.1 (| 2012 H is asso¬ 
ciated with the neutral torus and outflowing from the 
clouds as they are ablated. 

6. CONCLUSIONS 

We have described the rationale for defining a complete 
volume limited sample of bright local Seyferts, selected 
from the 14-195keV Swift-BAT catalogue. These AGN 
are complemented by a matched sample of inactive galax¬ 
ies, and we have shown that the host galaxy properties 
(stellar mass, morphological type, inclination, presence 
of a bar, distance) of these two samples exhibit similar 
distributions. Spatially and spectrally resolved observa¬ 
tions of the inner few hundred parsecs in these active and 
inactive galaxies will be used to study gas inflow and out¬ 
flow and the processes regulating it. The main points we 
have discussed here are: 

• A comparison of the properties of the Sy Is and 
Sy 2s shows that there is no bias for type in the 
selection. The only significant difference in AGN 
properties is the absorbing column Nh, the mean 
of which is two orders of magnitude greater for the 
Sy2s. This is consistent with unification schemes 
in which optical obscuration is greater for Sy 2s, 
but cannot distinguish between simple orientation 
dependence and instrinsic differences in torus prop¬ 
erties. 


• The fraction of Sy Is versus Sy 2s in this sam¬ 
ple, a nd also in the flux l imited 14-195 keV sam¬ 
ple of IWinter et all (1200911 . demonstrates that the 
true fraction of Sy2s at lo g Li 4-1 95 ~ 42. 5 -44 is 
50-60%, consistent with [Lawrence fc Elvisl (12010H 
and the spectroscopic subsample of lMerloni et all 
(I2014H . We show that the fraction of X-ray unab¬ 
sorbed Sy2s is at most a few percent and hence 
that, at these luminosities, optical obscuration and 
X-ray absorption usually occur together. 

• At higher luminosities log £ 44-195 > 44, while the 
optically obscured fraction drops rapidly, the X-ray 
absorbed fraction remains the same. We argue this 
implies that the inner boundary of the geometri¬ 
cally thick gas structure associated with the torus 
is roughly independent of AGN luminosity, having 
similar radial and height scales as for lower lumi¬ 
nosities. At low luminosities this gas structure is 
synonymous with the standard dusty torus; but at 
high luminosities it has an outer part which is the 
dusty torus, and an inner part which is a neutral 
dust-free torus that also contributes to the broad 
line region emission. In this context, the ‘receding 
torus’ model simply represents the location within 
this gas structure at which the predominant phase 
changes from ionised/neutral to molecular/dusty. 

• Finally, we note that the consistency of local (e.g. 
Swift-BAT) and more distant (e.g. the spectroscop- 
icall y classified AGN in the 0.3 < z < 3.5 sample 
of iMerloni et al.l I2014H samples suggests that the 
obscuring structure does not strongly evolve with 
redshift. 
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